Mass measurements of very neutron-deficient Mo and Tc isotopes and their impact on rp 

process nucleosynthesis 
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The masses of ten proton-rich nuclides, including the N=Z+\ nuclides ^^Mo and *^Tc, were measured 
with the Penning trap mass spectrometer SHIPTRAP. Compared to the Atomic Mass Evaluation 2003 a 
systematic shift of the mass surface by up to 1.6 MeV is observed causing significant abundance changes 
of the ashes of astrophysical X-ray bursts. Surprisingly low a separation energies for neutron-deficient 
Mo and Tc are found, making the formation of a ZrNb cycle in the rp process possible. Such a cycle 
would impose an upper temperature limit for the synthesis of elements beyond Nb in the rp process. 
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PACS numbers: 21.10.Dr, 07.75.+h, 26.30.Ca, 27.50.+e 

The astrophysical rapid proton capture process {rp pro- 
cess) [1] is a reaction sequence involving very neutron- 
deficient nuclides near the proton drip line, possibly up to 
mass number ~100. It is thought to cause astronomically 
observed X-ray bursts by repeated thermonuclear explo- 
sions in a thin, proton-rich fuel layer on the surface of a 
mass accreting neutron star. The nuclear energy generation 
in a burst and the composition of the burst ashes need to be 
reliably predicted to understand the properties of bursts and 
their impact on the structure of the neutron star crust, which 
is related to neutron star cooling [2]. The rp process may 
also produce the neutron-deficient ^-^^^^Mo and 96,98j^y 
topes, whose origin cannot be understood in the framework 
of standard nucleosynthesis [3]. Although it is not obvious 
that in X-ray bursts sufficient amounts of material escape 
the neutron star surface to contribute to galactic nucleosyn- 
thesis and though large rp process contributions may not 
be allowed by observed abundance patterns in meteorites 
10], such ejection is possible in principle [5]. Rapid proton 
captures on proton-rich nuclei are also an essential part of 
the up process occurring in the deepest, still ejected layers 
of core-collapse supemovae [6, 7], where interactions with 
neutrinos come into play as well. The process, however, 
depends sensitively on the astrophysical conditions and on 



the explosion mechanism. Therefore the reaction path is 
not uniquely defined, though calculations indicate that it 
can reach neutron-deficient isotopes in the Mo, Ru region 
and beyond. The contribution of these processes to the ori- 
gin of neutron-deficient isotopes in nature remains an open 
question. Improving the nuclear physics input is an impor- 
tant step to reliably calculate the isotopic production for a 
given astrophysical model. 

Nuclear masses strongly affect rapid proton capture pro- 
cesses ililil] as the low proton separation energies (S'p) en- 
countered near the proton drip line determine where photo- 
disintegration halts the sequence of proton captures and re- 
quires a /3+ decay, or a (n,p) reaction in case of the vp pro- 
cess, for the process to proceed. While much progress has 
been made in recent years with precision mass measure- 
ments of proton-rich nuclides using Penning traps, see e.g. 
|T3-l3l and references therein, the masses of crucial nu- 
clides for model calculations beyond Ak,1% have not been 
measured yet. Thus, one has to rely on the extrapolations 
of the Atomic Mass Evaluation 2003 (AME03) fl?] and on 
mass models. In this letter direct high-precision mass mea- 
surements of nuclides on the rp process path are presented. 

At the velocity filter SHIP [16] super- heavy elements as 
well as nuclides close to the proton drip line are produced 
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in fusion-evaporation reactions. For mass measurements, 
the exotic nuclides are transferred to the Penning trap mass 
spectrometer SHIPTRAP ifTvll . As will be demonstrated 
in this work, the combination of SHIPTRAP with the high 
beam intensities and superb separation capabilities of SHIP 
provides a unique reach for highly accurate mass measure- 
ments of extremely proton-rich nuclides. The products de- 
livered by SHIP are stopped in a helium-filled gas cell and 
extracted into an RF quadrupole. There, the ions are cooled 
and transferred in bunches to the Penning trap system. In 
the first Penning trap isobaric separation is carried out by 
mass-selective buffer-gas cooling, and in the second trap 
accurate mass measurements are performed using the time- 
of-flight ion-cyclotron-resonance technique (iB^. 

In the present experiment fusion-evaporation products in 
the element range from Rb to Tc were produced with a 
■^^Ar primary beam of 4 x 10^^ particles/second at ener- 
gies of 5.0 and 5.9 MeV/u reacting with a ^''Fe target with 
an area! weight of 0.45 mg/cm^, which has been enriched 
to 99%. The helium pressure in the gas cell was varied 
between 50 and 60 mbar. The mass measurements were 
performed using excitation times between 0.1 and 2.4 s. 
Mostly an excitation time of 1 .2 s was used, which corre- 
sponds to a mass resolving power (FWHM) in excess of 
10^ The masses of ten nuclides ^^Rb", »°^»^'»'^Sr, *^Zr, 
^^Nb, 85,86,87]y[Q ^jj^j s^Jq were measured with accuracies 
ranging from 4.6 x lO^'* to 2 x 10^'' using '^^Rb as refer- 
ence mass. Details of the experiment and the data analysis 
procedure will be the subject of a forthcoming publication. 

The results of the mass measurements are given in table 
U The masses of ^^Rb™ and 80,8i,84(gj. ^gj-g measured and 
found to be in agreement with the AME03. The more ex- 
otic nuclides, however, show a significant shift of the mass 
surface towards less bound nuclides as compared to the 
AME03 and thus follow the trend observed for Nb, Tc, Ru 
and Rh nuclides in previous experiments at SHIPTRAP and 
JYFLTRAP lilillj]. The masses of ^^Mo and *^Tc were 
measured for the first time and mass excesses are found 
to differ from the AME03 extrapolations by 1600 keV and 
1400 keV, respectively, well outside of the 300 keV ex- 
trapolation uncertainty given in AME03. Moreover, for 
86,87]y[Q^ which had only been indirectly measured before 
ifisll . the shift in the mass excess values exceeds the previ- 
ously given experimental uncertainties by several standard 
deviations. The masses of ^^Zr and ^^Nb also deviate from 
previous measurements listed in AME03 but agree with the 
values obtained at JYFLTRAP ifTlll . 

In order to investigate the influence of the new mass val- 
ues on astrophysical models, a consistent mass surface is 
required without artificial steps in the separation energies 
as would occur if only individual mass values were re- 
placed. For this reason, the new experimental mass data 
were added to the Atomic Mass Evaluation together with 
further recent experimental work [|Tol - IT3l . Taking into ac- 
count the entire set of experimental masses in this region, 
adjusted mass values were obtained following the proce- 
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FIG. 1. (Color online) Proton separation energies (Sp) of nu- 
clides with odd Z and even TV according to AME03 and AMEup. 
Filled symbols indicate measured values, open symbols indicate 
Sp where at least one of the two masses is extrapolated (#). 



dure employed in lfl9ll . For the nuclides measured in this 
work these values are given in table |l] column 6. Addi- 
tionally, a consistent mass surface requires updated mass 
extrapolations due to the large mass shifts found. Thus, 
a new local extrapolation in the region A=80-95 was per- 
formed using the method and programs of lll9ll . Together 
with the new set of evaluated experimental data and the 
previously reported AME03 extrapolated mass values for 
A<SO and yl>95 these data form an updated mass data 
set, which is labeled "AMEup" in this Letter. 

Even though the masses generally are shifted in the same 
direction, the AMEup table does also show changes in the 
separation energies compared to AME03. In particular, the 
Sp for ^'^Tc determined in this work experimentally for the 
first time is 868(6) keV, about 1000 keV lower than the 
previous extrapolated value of 1860(530) keV. Generally, 
one observes lower Sp for the updated mass surface and a 
more linear trend (Fig. [T}. 

To explore the impact of the new masses on the rp pro- 
cess in X-ray bursts, reaction network calculations using an 
X-ray burst model ||20|| were carried out. This single -zone 
model reproduces nucleosynthesis and energy generation 
reasonably well when compared to multi-zone burst mod- 
els but its computational speed allows one to explore nu- 
clear physics dependencies in detail. The baseline calcula- 
tion uses the nuclear masses of the AME03 and calculated 
Coulomb mass shifts ll2lll for nuclides beyond N=Z. The 
results are compared to network calculations based on the 
AMEup, combined with the same Coulomb mass shifts. 

The resulting final abundances show large differences 
between AME03 and AMEup in the region of y4=86-96. 
The largest change is found for A=86 where the abundance 
increases by a factor of 20 (Fig. |2]l due to the unexpected 
decrease in Sp of ^''Tc. The lower S'p of ^"^Tc changes 
the (7,p) and (p,7) rates such that the f3 decay branch for 
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TABLE I. Measured frequency ratios (col. 2) of the nuclides investigated (col. 1) with the corresponding mass excess values based on 
a ^'■^Rb mass of 84.911789738(12) u 1 15] (col. 3), as well as the corresponding AME03 values (col. 4), differences (col. 5), AMEup, 
(col. 6), Sp from AMEup (col. 7) and AME03 (col. 8) and their differences (col. 9). The symbol "#" indicates an extrapolated value. 
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The ground and isomeric states were not resolved. It is assumed that the isomeric state was populated, since it is favored by the reaction mechanism. 




60 70 80 90 100 110 



Mass number A 

FIG. 2. (Color online) Calculated overproduction factors (pro- 
duced abundance divided by solar system abundance) after an X- 
ray burst for the AME03 (dotted line) and the AMEup (solid line) 
mass sets. It has been assumed that all nuclei have decayed into 
the first stable isotope in each mass chain. The ratio of the over- 
production factor for a given isotope to the highest factor found 
for any isotope (10^) provides an upper limit for the contribution 
of X-ray bursts to the solar system abundance of that isotope. 



nuclides with mass number A=86 is strengthened. This 
change by a factor of 20 is by far the largest observed for 
abundances produced in rp process network calculations 
since the AME 2003 evaluation. It demonstrates that nu- 
clear physics uncertainties can be larger than estimated and 
can introduce large uncertainties in nucleosynthesis model 
calculations. Also the final abundance of A=9A increases 
by a factor of two. For mass numbers between 86 and 94 
the final abundances decrease. While the higher A=9A pro- 
duction is a step in the right direction, a much higher abun- 
dance would be required for A=9A to be co-produced suf- 
ficiently to explain the solar abundance of ^"^Mo (Fig. |2l). 
The change in A=9A production is due to the reduction of 
the Sp of ^^Ag (from 1040 keV to 790 keV), which is a 
result of the new mass extrapolation for ^^Ag. 

The results also demonstrate for the first time the exis- 



tence of a pronounced island of very low a-separation en- 
ergies (Sa) for neutron-deficient Mo isotopes, starting with 
^^Mo for which our measurement now provides an exper- 
imental Sa value (Fig. O. A similar situation is found for 
the Tc isotopes with our new Sa value for * ''Tc. This opens 
up the possibility for the formation of a ZrNb cycle induced 
by large ^^Mo(7,a) or ^■^Nb(p,a) reaction rates. The exis- 
tence of such a cycle had been discussed before 1230 based 
on the low Sa around '^^Mo predicted by the Finite Range 
Droplet mass Model (FRDM) |22]. However, the FROM is 
known to underpredict Sa severely in other isotopic chains, 
and the extrapolations provided by AME03 did not show 
such an effect. It was therefore concluded that the exis- 
tence of a ZrNb cycle is an artifact of the FRDM. Because 
of the new measurements this view has to be revised. 

The existence of a ZrNb cycle in the rp process is an 
important question. For a given Sa value of ^^Mo the cy- 
cle will form at a sufficiently high temperature, effectively 
providing an upper temperature limit for any rp process 
along the proton drip line to produce nuclei beyond A=84, 
including the light p-nuclei in the A=92-98 mass region. 
In order to explore whether our new results move this tem- 
perature limit into an interesting range, calculations in the 
Zr-Mo region were performed using a small test network 
with an initial *^Zr abundance. The degree of cycling was 
determined as a function of temperature and density. As 
a measure of the degree of cycling, the fraction of nuclei 
6cycio that end up in the A^=40 isotonic chain via ^*Moi'y,a) 
or ^■^Nb(p,a) and are not escaping the cycle via de- 
cay into A^=43 isotones were used. In order to explore 
the sensitivity of the results to nuclear masses additional 
mass data sets AME03x and AMEupx were created. Com- 
pared to AME03 and AMEup, a few masses of AME03x 
and AMEupx are changed by one standard deviation to 
the most favorable conditions for a cycle (low Sa and low 
Sp for ®'*Mo) to form. For each set of masses all forward 
and reverse reaction rates were recalculated using the sta- 
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FIG. 3. (Color online) (a) Alpha separation energies, Sa, for 
AMEup and the rp process path. Nuclides addressed in this work 
are indicated by a thicker frame, (b) Sa of molybdenum iso- 
topes for the AME03 (circles) and AMEup (squares). Filled sym- 
bols indicate measured values whereas open symbols indicate Sa 
where at least one of the two masses is extrapolated (#). 



tistical model code NON-SMOKER^'^^ version v5.8.1w 
123]. Compared to AME03, calculations with AMEup, 
AME03x, and AMEupx increase the ^'^Nb(p,a) reaction 
rate by factors of 40, 3 x 10^ and 8 x 10^ at 1 .7 GK, respec- 
tively, while leaving the *'^Nb(p,7) reaction rate largely un- 
changed. Nevertheless, even for AMEupx, the ^^Nb(p,7) 
reaction rate is still three orders of magnitude larger than 
^^Nb(p,a). Regardless of the masses, a cycle can there- 
fore only form when ^^Nb(p,7) is suppressed by a large 



inverse Mo(7,p) reaction, or when Mo(7,a) becomes 
significant. Only for AMEupx one finds a significant ftcycio 
for realistic temperatures up to 2 GK (Fig. |4). For higher 
temperatures, the radiation pressure in X-ray bursts would 
lead to expansion and cooling, and heavy nuclei such as 
^■^Mo would tend to be destroyed by photodisintegration. 
Both *'*Mo(7,a) and ^■^Nb(p,7) play a comparable role for 
the degree of cycling and the result does not depend sig- 
nificantly on density. While for high density more mate- 
rial is pushed past ^■^Nb, reducing the role of *^Nb(p,7), 
*'*Mo(7,a) turns out to be sufficiently fast to compensate, 
keeping the total ftcycic largely unchanged. A density range 
from 10^ to 5 X 10^ g/cm"^ was explored. 

Despite the significant cycling found for AMEupx 
masses at high temperatures, calculations with the full X- 
ray burst model, which reaches peak temperatures of 2 GK, 
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FIG. 4. Fraction of reaction flow branching into a ZrNb cycle for 
AMEupx (solid line) and AMEOBx (dashed line) as a function of 
temperature. AME03 and AMEup show no significant cycling. 



show that a cycle does not occur. The reason is that at the 
required high temperatures the reaction sequence already 
stops at ^^Ni. Temperatures that allow processing beyond 
^^Ni are only reached during burst cooling and are lower 
than required to form a ZrNb cycle. The formation is nev- 
ertheless a possibility in an environment where the temper- 
ature is rising slowly enough to enable the r-p process to 
proceed past ^^Ni before reaching high temperatures. An- 
other possibility would be an rp process with seed nuclei 
beyond ^^Ni. 

In summary, the masses of ten proton-rich nuclides were 
measured, among them for the first time the nuclides ^^Mo 
and ^"^Tc. These are the heaviest N=Z+l nuclides mea- 
sured so far, extending the previous measurements at other 
facilities. Thus this work marks a milestone towards high- 
precision mass measurements of the heavy N=Z nuclides, 
e.g. *"Zr and ^^Mo. The results show a shift of the mass 
surface towards less bound nuclei. The Sp of ^''Tc was 
experimentally determined for the first time and leads to 
large changes in the final abundance of ^=86 produced in 
the rp process. The abundance of A=9A, which is the pro- 
genitor for production of ^"^Mo, was found to be sensitive 
to the Sp of ^^Ag. First evidence is found of the existence 
of an island of low in the Mo region which is traversed 
by the rp process. This results in a revised, much lower 
prediction of the Sa value of ^"^Mo. Within uncertainties 
this opens up the possibility for the existence of a ZrNb 
cycle, which would impose an upper temperature limit for 
the rp process to synthesize nuclei beyond A=84. In order 
to quantify this temperature limit and to explore to which 
degree such a cycle poses a limitation to current models of 
the rp process precise mass measurements of ^°'®^Zr, ^'^Nb 
and ^''Mo are required. Both the low Sp of ^''Tc and the 
low Sa in the Mo region show that nuclear physics uncer- 
tainties in rp process calculations can be surprisingly large 
and that experiments are needed to put astrophysical mod- 
els on a solid foundation. 
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